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ADVANCED AIRFOIL DESIGN 
EMPIRICALLY BASED 

TRANSONIC AIRCRAFT-DRAG BUILDUP TECHNIQUE 

By W. D. Morrison, Jr. 
Lockheed-California Company 


SUMMARY 


Advances in airfoil section design applicable to aircraft optimized for 
transonic cruise offer improvements in the range factor M~, increases in 
wing thickness ratios , increases in wing aspect ratio , or reductions in 
wing sweep if applied to current transonic aircraft. 

To systematically investigate the potential of advanced airfoils in 
Advance Preliminary Design studies, empirical relationships have been derived, 
based on available wind tunnel test data, through which total drag can be 
determined recognizing all major aircraft geometric variables. This technique 
recognizes a single design lift coefficient and Mach number for each aircraft. 
Using this technique, drag polars can be derived for all Mach numbers up to 
+0.05 and lift coefficients -O.i^O to +0.20 from 

■^Design 

1 . INTRODUCTION 

In the preliminary design of advanced aircraft there is always the 
objective of incorporating into the design the most advanced technologies to 
evaluate their potential for future products . 

The concept of controlled supersonic flow development used in present 
day advanced airfoil design practices has evolved section characteristics 
e^diibiting improvements over more conventional airfoils for application to 
aircraft desigied for transonic cruise. The expression controlled supersonic 
flow development implies, that as the region of local supersonic flow develops 
and grows in extent over the airfoil chord, the shock wave terminating such 
a region remains weak. This phenomenon is in contrast to the very strong 
upper surface shock associated with more conventional design practices. 
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A mimber of prominent investigators have applied descriptive terms to 
their advanced airfoil design philosophy, i.e.. Dr. R. Whitcomb's super- 
critical airfoils, Pearcy's peaky airfoils, and Kom-Garabedian shockless 
airfoils. In addition, industry and the universities have developed their 
own advanced airfoil design practices and all indicate improvements over 
the conventional airfoils. All of these resulting airfoils may be considered 
supercritical; however, it has become customary to apply the term supercritical 
to those airfoils that in addition to having controlled supersonic flow develop- 
ment, aj-so carry a significant amount of aft loading due to an appreciable 
lower surface reflex near the trailing edge. Thus, of the above mentioned 
advanced airfoil designs, only those developed by Dr. Whitcomb of NASA fit in- 
to this more restrictive definition of the term supercritical. i 

At the present , limited nonsystematic 2- and 3-dimensional wind tunnel 
test data are available from these aforementioned sources. The explicit 
methods for the design and quantification of their characteristics reside with 
the principal investigators. Only the test results are available. 

This study has been undertaken, not in an attempt to devise an advanced 
airfoi'1 design procedure , but rather to collect available 2- and 3-dimensional 
test data under the generalized concept of "airfoils designed to a controlled 
supersonic flow development" and devise correlation techniques which will 
permit preliminary design evaluation of the potential of advanced airfoils 
for transonic flight applications. 

2. STUDY PREMISES AND OBJECTIVES 

State of the art and advanced airfoil 2-dimensional section shapes and 
approximate design condition press\ire coefficients are shown in Figure 1. The 
state-of-the-art airfoil shown here represents one of the sections tested during 
the development phase of the L-1011. The advanced airfoil shown is representa- 
tive of an outboard wing section employed on an advanced L-1011 wind tunnel 
model. Both sections are characterized by controlled upper surface supersonic 
flow development. However, the advanced airfoil rounder leading edge, flatter 
upper surface, and more reflexed trailing edge reduces the leading edge pressure 
coefficient peakiness, extends the near sonic flow further aft on the airfoil, 
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Figure 1. State of the Art and Advanced Airfoil Section Characteristics 






and produces a more highly aft loaded section; all contributing to a higher 
freestream drag divergence Mach number at a somewhat higher lift coefficient. 
Application of a similar pressure distribution to a 3-dimentional wing involves 
a design process which is complicated, undoubtedly varies between principal 
investigators , and to date employs empirical or other approximations 
(Reference l). In spite of these con^lications and probable variations, 
consideration was given prior to initiating this correlation in an attesgst to 
insure a degree of compatibility of definition and identification of parameters 
which might aid in the development or validation of theoretical programs. 

To this end, the assumption is made that a commonality exists between most 
theoretical design programs that a specified design-to pressure distribution, 
P/H, over the airfoil is stipulated at the critical or design lift coefficient 
and Mach number. This implies that for aJ.1 advanced airfoil design wind 
tunnel models there exists a somewhat common level of local Mach number over 
the airfoil at design conditions or for correlation purposes the same level of 
pressiire drag. The freestream Mach number at which these design conditions 
are met can be varied by the geometric parameters of sweep, aspect ratio, and 
thickness ratio. 

This suggests that the 3-dimensional wind tunnel test results can be 
correlated to 2-dimensional equivalence through relationships of design 
lift coefficient, design Mach n;miber, model geometry, and a relationship of 
the freestream Mach number to the design Mach number. 

3. APPROACH TO DATA CORRELATION 

Following the ass\imptions previously stated, and guided by the NACA work 
in developing data correlation techniques by means of the Transonic Similarity 
Rules (Reference 2), a data correlating procedure was evolved and is 
descriptively noted on Figures 2, 3, and i^. 

Two paths of data correlation are shown. Path I involves the deter- 
mination of design Mach number, design lift coefficient, and drag divergence 
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Figure 2. Path I of Data Correlation Technique - C . and M_ 

Derivation ^ ‘Design 
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Figure 3- Path I of Data Correlation Technique - Drag Divergence Mach 
Number Derivation 
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Figure 4. Path II Data Correlation Technique - Incremental Compressibility 
and Pressure Drag Derivation 




Mach niamber. Path II involves determination of the compressibility and 
pressure drag. 


Path I: Initial steps require the crossplotting of each wing-body 

data set as drag versus Mach number at constant lift coefficients; 
then computing M L/D for each data point and replotting as M L/D versus 
Mach number at constant lift coefficients. For maximum range, the 
lift coefficient and Mach n\mber at which the peak value occurs will, 
in general, be the airplane crviise or design condition. For practical 
reasons, i.e., flatness of the curve and to insure cruise at the 
highest speed, an arbitrary value of 0.99 (M L/D) was selected as 

the design point for M^ , and C,. 

Design . 

Design 


Using the same data crossplots of drag coefficient versus Mach number 
at constant lift coefficients, the three-dimensional drag divergence 
characteristics are then determined. To these data, AM corrections due 
to aspect ratio and sweep are applied to arrive at 2-D equivalent ' 
characteristics. Divergence has been defined as that Mach number at 
a constant lift coefficient at which the rate of change in drag 
coefficient with Mach number reaches a value of = 0.10. 

At the design C,. , Mach divergence, and design are synonymous. 

Xj 


• Path II: The basic drag polar, plotted as drag coefficient versus 

lift coefficient at constant Mach number, is selected for the deter- 
mination of the incremental compressibility and pressure drag con- 

tributions . Idealized induced drag is computed, as e = 1.0, 

and this quantity is removed from each polar. At the^owest test 
Mach number, at least a AM = 0.30 below design, the minimum pro- 
file drag is determined as the bucket of the curve and this drag 
quantity is removed from all polars. Referring back to Path I of 

the data correlation, . and C,. are selected and the 

Design 

Design 

adjusted polars are replotted at a constant AM from versus 

. The drag bucket of each of these polars is de- 

Design 

fined as the incremental compressibility drag and assumed to vary as 

a function of AM from M^ . , independent of lift coefficient. The 

Design 

remaining incremental drag above and below AC^ = 0, i.e., 

Design, 

is defined as the pressure drag and is assumed to vary as a function 

of AM from M^ . at a constant AC^ from C_ 

Design L 

Design 

The above da.ta reduction process was applied to each wing-body data set 
and the resulting data correlation is presented in the following sections. 
Where appropriate, 2-dimensional data were included. 


AC from C^. 

Xj Xj_ 
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h. DATA CORRELATION/EMPIRICAL RELATIONSHIPS 
U.l Data Base Summary 

The wing-body wind tunnel models incorporating advanced airfoil design 
sections are noted as through (see Figure 5); 2-dimensional data as 
through The L-1011 flight test and wind tunnel res\lLts , and W^, 

respectively, are included as a state of the art reference. is the 

L-1011 wind tunnel model test data incorporating an advanced airfoil section. 
Conventional wing design test data base, not noted on the figure but included 
in the correlation as another datum reference , was that obtained from 
References 3, and 5- throu^ test data were reported in References 
6 throvigh l8, respectively. The remaining references, 19 through 31, are 
included as a bibliography of advanced airfoil design tests. 

Wi , , and are the wing designs of Mr. L. R. Miranda of the Lockheed- 

California Company. embodies supercritical sections designed by Mr. J. A. 

Blackwell of the Lockheed-Georgia Company. Wg is the oblique wing of Dr. R. T. 
Jones of NASA-Ames with an airfoil section design by Dr. P. Garabedian. is 

a wing-body design by the Boeing Company of Seattle. The remaining wing and wing 
body models incorporate the supercritical wing sections of Dr. R. Whitcomb of 
NASA-Langley . 

The boxed-in numbers simply highlight the minimum and maximum geometric 
and design variables associated with this data base. 

k .2 Mach and C Drag Divergence 

Ij 

A comparison of conventional, state-of-the-art, and advanced airfoil 
design wing section drag divergence characteristics at a lift coefficient of 
0.50 is noted on Fig\ire 6. The data points represent either 2-D test results or 
3-D data corrected for the effects of aspect ratio and sweep using Table I. As 
an example, for an effective thickness ratio of 11 percent, advanced airfoil design 
gains over conventional airfoils are of the order of AM = +0.07- Compared to 
state-of-the-art airfoils, these gains are of the order of AM = +0.03 to 
+0.0it. Limiting thickness ratios, minimum and maximxmi, for incorporating 
advanced airfoil design practices have not been reported on to date and 
caution sho\iLd be exercised in extrapolating much beyond the available data. 
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TABLE I. INCREMENTAL MACH NUMBER CORRECTION TO DRAG DIVERGENCE 
DUE TO ASPECT RATIO AND SWEEP 


AR 

AM 

1.5 

0.097 

2.0 

o.o ?4 

2.5 

0.058 

3.0 

o.o 4 ? 

3.5 

o.o 4 o 

4.0 

0.035 

5.0 

0.028 

6.0 

0.023 

8.0 

0.018 

10.0 

0.014 

12.0 

0.012 

20.0 

0.008 



AM 

0 

0 

10 

0.004 

15 

0.008 

20 

0.016 

25 

0.028 

30 

0.042 

35 

0.057 

4 o 

0.075 

45 

0.095 

50 

0.115 

55 

0.137 

60 

0.160 
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Figure 5- State of the Art and Advanced Airfoil - Flight and Wind 
Tunnel Model Data Base Summary 
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Figirre 6, 2 D Approach - = 0.50 
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Recognizing that drag divergence Mach number may also he sensitive to 
other first order airfoil geometric variables, additional correlating para- 
meters were investigated and are noted on Figure 6. The parameter (M^ - l) 

2/3 ^ D 

versus (t/c) was suggested in Reference 2 and appears to improve the 

correlation. Effects of camber are also noted on this figure and reflect 
a slight variation in divergence Mach number with increased camber at this 
lift coefficient. At lower lift coefficients, the effects of high camber 
become more predominant. For purposes of this study, a linear variation of 
- l) vs (t/c) , representative of cambers up to a;2.2^, is assumed 

and this correlation is presented for lift coefficients from 0.10 to 0.80 on 
Figures 7 through lU. Tables II and III are included to aid in determination 
of the Mj^ 2 relationships to the chosen parameters (1^2 d ~ 

(t/c)^/^. 

H .3 Design Lift Coefficient 

The determination of Ct„ . is the key to this data correlation technique. 

•‘■'Design 

Through this definition there is a single optimum polar selected for a given 
aircraft geometry. 

The initial study premise, that each wing-body wind tunnel data set has 

its unique design lift coefficient, implied that certain configuration 

relationships would determine' this lift coefficient. The classical definition 

of Ct , , (the lift coefficient corresponding to (L/D ) suggested those 
■‘■'optimum ^ max 

correlating parameters, since ^Lopt ~ ^ ( V eTTAR \/c^. ) On Figure 15 the 
CLDesign wing-body configurations of this study are noted as a function 

of wing aspect ratio. By normalizing those data to a single level of = 
0.0160 and recognizing section camber as a contributor to "the 

relationship shown on the left of Figure 15 was determined. 


h.U Compressibility and Pressure Drag 

An assessment of the magnitude of the various incremental contributors 
to total drag is in order prior to review of the following section. 

The percentage contribution of the compressibility plus pressure drag 
to the total drag of an existing wide-body transport at the design conditions 
of M = 0.85 and 0.46 lift coefficient is some l 4 percent. This contribution 
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Figure 8. Correlation - = 0.20 




Figure 9 . M^„vst/c Correlation - =0.30 
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Figure 11. 2 p t/c Correlation " = 0.50 
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Figure 12. M vs t/c Correlation - C = 0.60 
D C--D L 


- 0.7 
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vs t/c Correlation - C_ = O.70 
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Figure l4. ^ 2 D Correlation - = 0.80 


table II* 


^ 2-B 



X) 2-D’ 


and 






% 2-D 

0.60 

0.62 

0.6U- 

0.66 

0.68 

0.70 

0.72 

0.7^ 

0.76 

0.78 

0.80 

0.82 

0.84 

0.86 

0.88 

0.90 


2 

% 2-D 

0.3600 

0.3844 

0.4096 

0.4356 

0.4624 

0.4900 

0.518^ 

0.5476 

0.5776 

0.6084 

0.6400 

0.6724 

0.7056 

0.7396 

0.7744 

0.8100 



-0.6400 

-0.6156 

-0.5904 

-0.5644 

-0.5371 

-0.5100 

.0.4816 

-0.4524 

-0.4224 

-0.3918 

-0.3600 

-0.3276 

-0.29^^ 

_0.26o4 

-0.2256 

-0.1900 
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TABLE III. 


t/c, (t/c)^/^, (t/c)^/^, and (t/c)^/^ RELATIONSHIPS 


t/c 

(t/c)5/3 

(t/c)^/3 

(t/c)^/3 

o.o 4 

0.00468 

0.1170. 

0.3420 

0.05 

0.00679 

0.1357 

0.3684 

o.o 6 

0.00920 

0.1533 

0.3915 

0.07 

0.01189 

0.1699 

0.4121 

0.08 

0.01485 

0.1857 

0.4309 

0.09 

0.01807 

0.2008 

0.4481 

0.10 

0.02155 

0.2154 

0.4642 

0.11 

0.02526 

0.2295 

0.4791 

0.12 

0.02919 

0.2433 

0.4932 

0.13 

0.03337 

0.2566 

0.5066 

o.i 4 

0.03775 

0.2696 

0.5193 

0.15 

0.04235 

0.2823 

0.5313 

o.l 6 

0.04716 

0.2947 

0.5429 . 

0.17 

0.05218 

■0.3068 

0.5540 

0 .l 8 

0.05739 

0.3187 

0.5646 

0.19 

0.06280 

0.3304 

0.5749 

0.20 

0.06841 

0.3419 

0.5848 

















is essentially equally divided as 7 percent compressibility and 7 percent 

pressure. The remaining drag is comprised of 33 percent induced and 53 

percent friction, dirtyness , leakage, and interference. As the Mach number 

is increased to 0.89, AM = +0.0U, the {AO + AC_. ) contribution increases 

JDc up 

to some 30 percent. Reducing the Mach niomber to M = 0.80, AM = 0.05, this 
contribution is in the order of 10 percent of total drag. 

In evaluating the compressibility plus pressure drag of 2-dimensional or 
3-dimensional wing-body test data, consideration must be given to the test 
Reynolds number, location of the boundary layer trip grit, and wind tunnel 
blockage. The sensitivity of + ACj^^to Reynolds number is noted on 

Figure l6. The design conditions were approximately Cjj = 0.60 at a Mach number 
of 0.79- A.S the test RN is increased, the level of ( A Cj)^ + Cj)^) tends to 
reduce or bucket as the design Mach number is approached. These effects will 
be noted in the following 3-D correlations and are in the main evidenced in the 
ACj)^ term. 

U.U.l Compressibility Drag. - On Figure 17 is presented as a function 

of AM from 'the wing-body data sets and the L-1011 state-of-the-art 

reference flight test and wind tunnel data. At the top of the figure , the 

unadjusted compressibility drag shows scatter on the order of 20 counts of 

5/3 

drag, i.e., “ 0.0020. The correlating factors of (t/c)^ and camber 

h/c % (lower curve, Figure 17) were those suggested in NACA TR 1253 - Cor- 
relation by Transonic Similarity Rules. Up to ^ggj^gj^ correlation of data 
appears exceptionally good. Data scatter above ^ggj^gj^ attributable to 

varying test Reynolds number, test grit location, camber, and section variation 
with span, or to the degree of refinement of model following flow visualization 
studies of shock buildup across the span. A fairing of these data against a 
logorithmic scale is presented on Figure l8. 

h.k.2 Pressure Drag. - Pressure drag is presented on Figure 19 through 28 
for AC 's ranging from minus 0.40 to plus 0.20. Figure 2h (upper c\irve) 
presents the uncorrelated data sets at ^Lpgg^g^- Design C^ for the data sets 
ranged from 0.275 to 0.55 and Mach design from 0.73 to 0.97. Approximately 
20 counts of data scatter are noted up to design Mach nijmber with increasing 
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(9'0=IAI ® Pd-Pd) 



Figure l6. 2-D Compressibility and Pressure Drag Variation 

with Reynolds Niamber 
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■ 0.20 - 0.16 - 0.12 - 0.08 - 0.04 0 0.04 

AM=Moo-MdeS. 

.Figure I8. Fairing of Incremental Compressibility Drag - Summary 
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Figure 19 . Incremental Pressure Drag Variation with 
Incremental Mach Number from Mach Design 

Ap = -0.40 
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AM = Moo - Mpgg 

Figure 20. Incremental Pressure Drag Variation with Incremental Mach Number 

from Mach Design -AC_ = -O. 3 O 
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Figui"® 21. 


AM=Moo -MdES 


incremental Pressure Drag 

from Mach Design 


Variation 

- 0.20 


with 


incremental Mach Number 
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(1 +i:%/io) 



AM = IVloo 


Figure 22. Incremental Pressure Drag Variation with Incremental Mach Number 

from Mach Design -AC = -0.10 
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(1 +il%/10) 



“ Moo * 

Figure 23. Incremental Pressure Drag Variation with Incremental Mach Number 

from Mach Design -ACt = -0.05 

L 


33 



AM = Moo - Mpgg 

Figure 24. Incremental Pressure Drag Variation with Incremental Mach Number 

from Mach Design -AC^ = 0 
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AC 



AM — Moo * 

Figure 25. Incremental Pressure Drag Variation with Incremental Mach Number 

from Mach Design -Ac = +0.05 
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(1 +r%/10) 



- 0.20 - 0.16 - 0.12 - 0.08 - 0.04 0 0.04 0.08 

— Moo - 

Figure 26. Incremental Pressure Drag Variation with Incremental Mach Number 

from Mach Design -ACt = ^-0.10 

L> 
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(Ol/% 1) 



“ Moo ’ ^QES 

Figure 27. Incremental Pressure Drag Variation with Incremental Mach Number 
from Mach Design - An = +0.15 

Jj 
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Q. 



< 




AM = Moo -Mpgg 


Figure 28. Incremental Pressure Drag Variation with Incremental Mach Number 
from Mach Design -Ac_ = +0.20 

Lf 
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scatter above. Again, the correlating factors siaggested in NACA TR 1253, i.e., 
1/3 

(t/c) and camber h/c were en 5 >loyed to adjust the data sets for varying 
geometries (lower curve. Figure 2h) . Again the correlation is improved up 
to design Mach number with increasing data scatter above. 

Around . , grouping of the data pinpoints the desirable reduction in 

pressure drag evidenced from the 2-D tests noted on Figure l6 . Figures 19 
through 21 and 25 through 28 present both \mcorrected and adjusted ACp 

P 

for AC^'s ranging from -O.UO to +0.20. (Note scale changes.) At a 

AC_ = -O.UO, Figure 19, and AM's above NL , the effects of camber do not 
L ° Design 

appear to be sufficiently accomted for by the h/c correction. The deviation 

of the data may be attributable to lower surface shocks and separation which 

was also evidenced in the Mach divergence characteristics of Figure 7 through 

lU . It should be pointed out, however, the a AC = -0 .kO represents a very 

low lift coefficient, i.e., 0 to +0.15, for the majority of the data sets and 

is seldom a lift coefficient corresponding to any flight phase of a transport. 

Figure 29 summarizes the ACn faired variation with k M for AC , of -0 .kO 

p Jj s 

to +0.20. 


U.5 Zero Lift Pitching Moments 

Large negative Cj^^ has been associated with advanced airfoil design 
practices and raises a concern that this might impose significant trim 
penalties . 

On Figure 30 the variation of Cj,^ with AM indicates a fairly constant 

value up to ^-ith some reduction above. The upper curve of Figure 30 

indicates at M^ . that C™ correlates" as might be expected with camber. 
Design ™o 

For preliminary design purposes, the assumption may be made that under 

trimmed'^conditions and for configurations incorporating aft tails , the wing 

must carry an additional lift equal to but opposite in sign to the tail down 

it 

load. This increment in lift is AC^ = -Crr, /~=-. At the design C^ and Mach 

^ ° ^ tt ^ 

number for a 10^ thick airfoil cambered 1.6^ with an — =• = 3.0, the wing will 

c 

carry approximately 3 counts of trim drag to provide extra AC . The increased 

Li 

down wash at the tail and resulting angle of the tail lift vector may subtract 
from this trim drag; however, this increment would be too difficult to predict 
in the normal preliminary design study. 
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Figure 29- Fairing of Incremental Pressure Drag - Summary 
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Figure 30. Zero Lift Pitching Moment Variation with Camber 
and Incremental Mach Number from Mach Design 


iil 


h.6 Buffet Onset 


For purposes of this discussion, buffet onset is defined as the break in 
axial force versus lift coefficient or for those tunnel models employing root 
strain gages as the lift coefficient corresponding to a rapid rise in the RMS 
root bending moments . 

A systematic buffet investigation was conducted by NASA and results were 

reported in NASA -TN 5805. The wings of this test series employed 63, 6k, 65, 

i+5°, and aspect ratios U. to 6. Previous buffet correlation attempts by 

Lockheed employing the data of this report had shown a reasonably good 

correlation of buffet onset with Ct^ . . Some buffet data were obtained on 

^Design 

the advanced airfoil designs employed in this study, and a comparison was made 
(Figure 31) to compare buffet characteristics of conventional and advanced 
airfoil sections. 

The variation of the advanced W and W wing/body C Buffet versus AM 

( XU 1j 

exhibits a vastly different character. exhibits the same variation with 

AM as a conventional airfoil, while W^, which is more representative of 
actual flight experience, has a very different variation. The reason for 
these differences in buffet onset variation are not understood and certainly 
additional effort is required in buffet prediction techniques. 


1+.7 Airfoil Section - Pressure Drag - Form Factor 


The friction drag buildup process for preliminary design configurations 
will generally involve knowledge of both the wing and body form factors as a 
function of wing section and body geometry. 



F.F 


S 

S 


Wet 

Ref 


On Figure 32 this form factor has been developed for advanced airfoil 
sections and also presents data for a NACA 65 series section for reference. 
The parameter F.F = Cp^/2 is simply the 2-dimensional minimum section drag 
coefficient divided by twice the flap plate skin friction coefficient at the 
test Reynolds number. The factor of 2 accounts for both upper and lower sur- 
faces. The data, from top to bottom on Figure 32, is from 2-D tests on a 
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Figure 31- Lift Coefficient for Buffet Onset 
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Figure 32. Conventional State of the Art, and Advanced Airfoil Section - 
Ratio of Minimum Drag to Theoretical Skin Friction Drag 
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conventional NACA 65 , 213 a = O .5 airfoil, ^ percent thickness ratio state-of- 
the-art airfoil, and 10 percent and 11 percent thickness ratio advanced air- 
foils, respectively. The flagged versus unflagged symbols represent same 
model tested in two different facilities. On Figure 33 j the average fairing 

of Cj) /2 C„, is shown as the form factor versus the section thickness ratio. 

o /■ 

The NACA 65 series airfoil results confirm the variation given in the HAS data 
sheets and USAF DATCOM. At a thickness ratio of 10 percent the advanced air- 
foil appears to carry an approximate 10 percent increased subsonic pressure 
drag over the conventional airfoil sections. Subsequent advanced airfoil test 
results (data unpublished) of 2nd generation designs are included bn Figure 33* 
Additional information on 2nd generation airfoils is included in Appendix A. 

Figures 3^» 35) and 36 , which were not derived as part of this study, 
are presented as part of the total drag buildup procedure to be discussed in 
a later section. 

5 . drag buildup procedures 

The correlation techniques developed in the previous sections can now 
be combined to define a systematic procedure for the total drag buildup of 
arbitrary transonic aircraft incorporating advanced airfoil wing sections. 

The geometric characteristics of an example twin-engine transonic trans- 
port are itemized on Figure 37 and will be used to demonstrate the drag build- 
up process. It should be emphasized here that, any interference or compress- 
ibility dra^ of nacelles or tails, and incremental drag due to inlets, antennas, 
or other proturbences, is an additive drag element that must be derived from 
other sources. 

The initial step in the drag buildup procedure, graphically described on 

Figure 38 ) involves the determination of the aircraft design lift coefficient 

and design Mach number. The derivation of *"Lpesign determination 

of the minimum pressure drag at M = G .60 at the nominal mission altitude . For 

this example 30,000 ft has been selected. Aircraft component drag buildup is 

on Figure 38 along with those text figures required for the computation. Using 

this computed Cp. . , percent wing camber, and wing aspect ratio; Cj is 

^pmin -Resign . 
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, Figure 36. Compressibility Correction to Skin Function Coefficient 


AIRCRAFT GEOMETRY 



Figure 37 - Example Twin Engine - Advanced Transonic Transport Geometry 
















determined. Computation of Mach design involves determination of the 

2-dimensional equivalent divergence Mach number for the wing section at 

the design lift coefficient. This M^^ ^ then corrected for wing aspect 

ratio and sweep to obtain the 3-diniensional drag divergence Mach number, which, 

at Ct , is equivalent to 

^Design Design 

At this point, the lift coefficient and Mach number have been determined 
for best cruise. The total drag polars can now be defined for any Mach- 
Altitude condition; however, for purposes of this example buildup, the best 
cruise Mach number is chosen. Note that if other than best cruise altitude 
is selected for off design polar determination, minimum pressure drag must be 
coB 5 )uted based on the Reynolds number and associated with that altitude 
Mach number condition. In the final aircraft preliminary design sizing proc- 
ess, off design Mach number polars would also be investigated to assure an 
airframe-propulsion system match for all mission segments. 


The next step in the drag buildup procedure involves the computation of 

the minimum pressure drag, the induced drag, compressibility drag, pressure 

drag, and trim drag at the design Mach number of O .85 and altitude of 30>000 

ft. These computations are shown on Figure 38 in tabular form. Most of the 

steps are self evident. Coliimn is equivalent to the trimmed 0,. = W/ 

Jj KJii? 

The AC due to trim is an additional lift which must be carried on the wing 
Jj 

to balance the tail down load required to trim the wing pitching moment. For 
this example a zero static margin was selected; hence, the required trim is 
only that amount necessary to overcome the negative c^ of the wing = -O.OU. 
Perhaps, with advancement in active control technology, this trim requirement 
on resulting drag can be reduced to zero at the design conditions. Column © 
is the incremental lift that the wing is experiencing from the design lift 
coefficient. The remaining computations are relatively straightforward, with 


the possible exception of ( 17 ; and (20j . In determining the AC^^ , the param- 
eter presented on Figure 29 must be crossplotted versus AC^^ at ^AM = 0 to 
allow interpolation of the pressure drag contribution at the AC, represented 

ij 

by column © . The final drag represented by the summation performed for 
column must be plotted versus the C, = W/ of column ©to represent the 


L ' qs 

trimmed aircraft polar. For determination of drag polars at other Mach-Altitude 
conditions, the same process is followed, i.e., con^jutation of columns © 
through j respectively. 
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6. 2-D AND 3-D COMPARISON OF COMPRESSIBILITY AND PRESSURE DRAG 

The objective of this section is to establish an equivalence comparison 
between 2- and 3-dimensional compressibility and pressure drag which might 
permit (l) an assessment of an ultimate potential of advanced airfoil design 
practices, (2) indicate some of the- penalties associated with applying 2-D 
characteristics to the 3-D wing, and (3) indicate a degree of data sensitivity 
to some of the major test variables . 

On Figure 39 the fairing of the 3-d.imensional wing-body compressibility 
and pressure drag increments; adjusted for the effects of aspect ratio, sweep, 
camber, thickness ratio, *^LDesign design’ the methods developed in 

this report; are compared with 2-dimensional test data for a 10 percent thick- 
ness ratio, 1.64 percent camber, advanced airfoil section. The flexibility of 
the definition of advanced airfoils; sensitivity of the component drag data to 
the test variables; varying design techniques employed by the principal inves- 
tigators in arriving at the 3 -dimensional wing section, camber, and twist; and 
the varying degree of model refinement during testing all contribute to the 
data scatter previously evidenced in the 3-dimensional data. However, the 
conservatism in the faired 3-D data relative to the 2-D data appears mainly 
attributable to differences in test Reynolds Number, i.e., 3-D wing-body tests 
are generally conducted at RN = 2 to 3 million. Additional coit^iarisons are 
noted on Figure 40 for other AC , . The agreement between 2- and 3-D becomes 
increasingly better at the lower lift coefficient with increasing disagreement 
above Design C^. 

At the present stage of development of advanced airfoil design practices, 
it is believed that for preliminary design of aircraft for transonic cruise 
applications the potential of advanced airfoils is best represented by the 
3-dimensional data fairings of the previous sections. Additional tests and 
data analysis, the subject of the next section, certainly are in order to 
sort out the 2-dimensional and 3-dimensional equivalence differences and in 
particular the sensitivity to Reynolds Number. 
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7 . REC OMMEUDATIONS 


The advanced airfoil data correlation and drag buildup technique arrived 
at in this reporting may assist in determining the potential of advanced air- 
foil^ applications ; however, there remain major and significant areas of fur- 
ther study. These study areas are itemized below: 

1. Systematic definition of 2-D design-to criteria and testing to 
isolate the sensitivity and limitations of the geometric varia- 
tion on 

• design lift coefficient and Mach number 

• lift coefficient for buffet onset 

• maximum lift coefficient 

• zero lift pitching moment 

• minimum pressure drag 

2. Determination of the 2- and 3-dimensional equivalence in the above 
parameters and 3-D design procedures using 2-D data. 

3. Systematic 3-D testing of higher aspect ratio wings including buffet 
onset studies. 

4. Continued development of empirical relationships such as those devel- 
oped in this study extending into maximum lift and pitching moment 
correlation . 

5. Empirical correlation of lift, drag and moment characteristics of 
high- lift devices on wing/body's employing advanced airfoil 
sections. 

6. Systematic wind tunnel investigations to determine the potential of 
advanced airfoil design practices for supersonic cruise applications. 
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APPENDIX A 


Subsequent to the development of the drag buildup procedures presented 
in the body of this report, additional advanced airfoil wing-body wind tunnel 
data was made available through NASA Langley. (See Figure A-1. ) 

Analysis of these data by the previously discussed procedures resulted 
in appreciable reduction in the parameters and ^^j)p some slight 

penalty in drag divergence Mach number. 

These new fairings of qj ^^Dp’ '^^Dc noted on Figures A-2 

through A-7 and are recommended for use as "Advanced Airfoil Design-Potential" 
drag buildup. 

The more salient features of these advanced airfoils are 

• A reduction in the compressibility drag parameter, i.e., drag creep 
for all Mach numbers. The fairing of the data on Figure A-3. reflects 
levels attained on 2-dimensional sections, i.e., 

• The pressure drag parameter reflects appreciable reduction in pressure 
drag and the Mach number is increased towards Mach design. (See 
Figures A-3 through A-7. ) 

• Drag divergence Mach numbers are reduced approximately AM = -.015 
over the Advanced Airfoil reported on in the body of this report. 

(See Figure A-2.) 

• The form drag appears to be unchanged. Figure 33 form factors are 
still recommended. 

• For a wing-body having an effective wing thickness ratio of 10 per- 
cent, the "Advanced Airfoil Design-Potential" compressibility and 
pressure drag would be approximately l4 drag count lower at 

and C]j)esign Advanced Airfoil reported on in the body of 

this report. 


A-1 




Q) 

W 

c6 

m 

aj 

-P 

c6 

Q 


0 

m 

1 

bO 

C 

•H 

& 



A-2 



D2-D 
















